Laser micro-fabrication of concave, low-roughness features in silica 



D. Hunger, 1 C. Deutsch, 2 R. J. Barbour, 3 R. J. Warburton, 4 and J. Reichel 2 '* 

1 Max- Planck- Institut fur Quantenoptik and Fakultat fur Physik der 
Ludwig- Maximilians Universitat, Schellingstrafie 4, D 80799 Miinchen, Germany 
2 Laboratoire Kastler Brossel, ENS/UPMC-Paris 6/CNRS, 24 rue Lhomond, F-75005 Paris, France 
3 School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, UK 
^Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland 

(Dated: September 26, 2011) 

We describe a micro-fabrication method to create concave features with ultra-low surface rough- 
ness in silica, either on the end facets of optical fibers or on flat substrates. The machining uses a 
single focused CO 2 laser pulse. Parameters are chosen such that material is removed by thermal 
evaporation while simultaneously producing excellent surface quality by surface tension-induced 
movement in a low- viscosity melt layer. A surface roughness a ~ 0.2 nm is regularly obtained. The 
concave depressions are near-spherical close to the center with radii of curvature between 20 and 
2000 /im. The method allows the fabrication of low-scatter micro-optical devices such as mirror 
substrates for high- finesse cavities or negative lenses on the tip of optical fibers, extending the range 
of micro-optical components. 



CO2 laser light is a powerful tool for the surface treat- 
ment of fused silica. Strong absorption of the 10.6 fim 
light within the first few jim enables controlled surface 
melting using moderate intensities. A central benefit of 
such treatment is that surface tension in the molten layer 
smoothens out surface roughness up to a scale compara- 
ble to the thickness of the melt layer. This was recog- 
nized as an efficient way to polish optical surfaces [1, 2]. 
Owing to the surface- minimizing effect of the surface ten- 
sion, melting of large-scale domains which penetrate con- 
siderably into the volume will result in the formation of 
convex structures. While this behaviour has been suc- 
cessfully employed for the fabrication of convex micro- 
optical structures [3-6], it is undesired in the context of 
surface polishing. At higher intensities, material removal 
by evaporation becomes important. In the context of op- 
tical surfaces, evaporation has mainly been studied as an 
unwanted side-effect of damage repair [7, 8] and polishing 
[2]. Nevertheless, it holds promise for micro- machining 
of optical components [9]. 

Here we describe micromachining of concave profiles 
on silica surfaces in a simple process using a single laser 
pulse. We work in a regime where surface evaporation is 
dominant and shapes a depression, while melting is re- 
stricted to a thin layer that smoothens the surface on a 
short scale. This avoids the need for scanning techniques, 
special atmospheres, preheating, or polishing pulses. The 
method allows fabrication of microscopic concave sur- 
faces on the tip of optical fibers as well as on bulk fused 
silica substrates. It was originally developed to realize 
high-finesse fiber Fabry-Perot cavities (FFPs) that en- 
abled the first cavity quantum electrodynamics (CQED) 
experiments with Bose-Einstein condensates [ ], and 
very recently has been used with single organic molecules 
[11] and quantum dots [ ]. More details on these FFPs 
can be found in [13, 14]. Here, we present an investiga- 
tion of the physical process leading to the formation of 
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FIG. 1: (a) Surface profile of a machined fiber measured by 
the optical profilometer. Beam parameters are P = 420 mW, 
r = 27 ms and w = 28 /xm. (b) Cut through the center of the 
profile shown in (a) (red solid line) and its fit to a Gaussian 
(black dash-dotted line), (c) Central part of the same data 
(red solid line) along with a circle fitted to the center yielding 
R — 70.8 /im (black dash-dotted line). Also shown is the local 
radius of curvature as calculated from a high-order polynomial 
fit to the data (green dotted line). 



these ultrasmooth concave depressions. 

We use an RF-pumped CO2 laser (Synrad Firestar 
v20) with a nominal power of 20 W. Pumping is pulsed 
at 20 kHz, the average power P being controlled via the 
duty cycle. Laser power fluctuations are minimized by 
reducing unwanted back-reflections to the laser with a 
quater-wave plate and polarizer, and we obtain energy 
fluctuations of 1.2% for a 4 ms pulse train. The beam is 
focussed onto the sample surface. To position the sam- 
ple with respect to the focal waist, we use a long working 
distance microscope that allows us to inspect the surface 
via a dichroic mirror. We estimate the alignment preci- 
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sion to 2 /im transversally and 7 /im along the beam axis, 
limited by the depth of field of the microscope. 

Our substrates are metal-coated, single-mode or 
gradient-index optical fibers (Oxford Electronics) with 
125 and 200 /im cladding diameter, and 2 mm thick fused 
silica glass plates (Hellma QS). The fibers were not 
cleaned after cleaving; the plates were roughly cleaned 
with ethanol. After machining, surfaces were analyzed 
with several methods. Surface roughness <r rms is evalu- 
ated by atomic force microscopy (Dimension 3100 AFM) 
and optical loss measurements, discussed further below. 
The overall shape is measured by optical interferometric 
profiling (Fogale Micromap 3D). Figure 1 shows a typical 
machined surface. We observe profiles with an approxi- 
mately Gaussian shape. The geometrical properties used 
for characterization are the central radius of cuvature i?, 
the total depth t and the structure diameter d. We eval- 
uate R from a polynomial fit to the center of the struc- 
ture. We define t as the height difference between the 
maximum and the minimum elevation, and d as the dis- 
tance between the turning points from convex to concave 
curvature of a high-order polynomial fit. 

We machined a large set (several hundreds) of fibers, 
and a series of glass plates, with a broad range of pa- 
rameters. Laser powers were between P = 300 mW and 
2 W at the sample surface and pulse trains had durations 
r = 4 to 120ms (consisting of 80 to 2400 pulses). Beam 
waists were between w = 21 and 93 jam. With these pa- 
rameters, the resulting structures are t = 0.01 ... 4 /im 
deep and have diameters of d = 10 ... 60 /im, leading 
to radii of curvature R = 20 . . . 2000 /im. Fig. 2 shows 
how the results depend on the laser parameters. For in- 
creasing r, the structures become deeper and R becomes 
smaller. For increasing w, both d and R increase also. 
The measured R agrees well with the value R g = d 2 / (St) 
expected for a gaussian shape (inset in Fig. 2): of the 
three parameters {d, t, R}, only two are independent and 
fully determine the third. 

The fiber data show a spread for fixed parameters of up 
to 13% in the profile depth, 23% in the profile diameter 
and 13% in the radius of curvature. Probable causes for 
this spread are the fluctuations in the laser energy and 
the limited positioning precision along the optical axis. 

The possible range of laser parameters for the fabrica- 
tion of concave profiles is restricted by two physical lim- 
its. For very short pulses at high intensity, melting and 
resolidification occur faster than the surface smoothing 
process [ ]. (Our parameters do not reach this limit.) For 
very long pulses, the melt layer extends far into the vol- 
ume and leads to global contraction into a convex shape. 
This regime is announced by the formation of a bulge 
around the depression at our longest pulse durations (red 
dashed curve in Fig. 3 (a)). The bulged depressions de- 
viate from a Gaussian shape, becoming shallower and 
broader. 

Some insight into the evaporation process can be 
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FIG. 2: Experimentally realized geometries: R as a func- 
tion of t. Beam powers and waists: ■ 600 mW, 26 /im. 
□ 1.85 W, 73/xm. • 373 mW, 27 /mi. o 0.9W, 63 /mi. 
A 540 mW, 28 /mi. A 575 mW, 63 /mi. Data for a given pa- 
rameter set were taken by varying r; r increases from left to 
right. Squares: bulk material, circles: 0125 /im fiber, trian- 
gles: 0200 /im fiber. Inset: The measured R agrees well with 
the value R g calculated from d and t. (Straight line: R = R g .) 
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FIG. 3: (a) Profile as a function of r. Shown are cuts through 
the center of 0125 /im fibers for P — 572 mW, w — 51 /im. 
The curves are vertically offset by 100 nm for readability, (b) 
Central depth and (c) diameter of the depression as a func- 
tion of r for similar beam parameters on different targets: 
• 0125 /im fibers. A 0200 /im fibers. ■ 2 mm thick fused 
silica plate. Solid lines: Fit to the described model, where 
k is used as fitting parameter to the data in (b). We find 
k 2.05 W/mK for the 125 /im fibers, k 2.0W/mK for the 
200 /im fiber, and k = 2.55 W/mK for the glass plate. The 
same parameters are used for the model curves in (c) , yielding 
a reasonable agreement for the plate, but not for the fibers. 
Beam parameters: • 600 mW, 26 /im. A 540 mW, 28 /im. 
■ 600 mW, 26 /im. 



gained from a simple model in which the surface tem- 
perature profile is used to calculate the ablation rate. 
The local temperature is calculated considering the 
absorption of a laser beam (intensity profile I(r) = 
Iq exp(— 2r 2 /w 2 )) by the surface of a half-space filled with 
a homogeneous medium. The medium is characterized by 
its absorption coefficient A, thermal conductivity ft, and 
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thermal diffusivity D. (Numerical values for fused silica 
at room temperature are A = 0.85, k = 1.38 W/mK, and 
D = 7.5 x 10 _9 m 2 /s.) The resulting surface temper- 
ature profile is obtained from the two-dimensional heat 
flow equation [7, 15]: 



T(r,r) = 
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%vd r 
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dr'. (1) 



The temperature in the center can be evaluated ana- 
lytically, T(0,r) = AI w/(V2ttk) arctan y^r/(w 2 /SD). 
The material removal rate (velocity of the evaporation 
front) is related to the local temperature as v(r,r) = 

Voe -U/k B T(r,r)^ where u = 3 6e y ig the latent heat of 

evaporation per atom and = 3.8 x 10 5 cm/s [ ]. It be- 
comes significant above ~ 2000° C and has a very strong 
temperature dependence in this temperature range. This 
simple expression neglects all non-linearities as well as 
the boundary condition set by the finite size of the fiber 
that modifies the transverse heat conduction. 

For short pulses (no significant bulge), the depth of the 
structures is reproduced well by this model (Fig. 3(b)). 
We treat n as a free fitting parameter on account of its 
strong temperature dependence; the resulting k values 
are in the range of published values [16] for hot silica. 
(The small difference between the k values for the two 
fibers is not significant considering the uncertainties of 
the power and waist measurements.) The observed pro- 
file diameters for plate machining are in reasonable agree- 
ment with the model. For fiber machining, however, the 
profile diameters are systematically larger than predicted 
by the model and also show a different parameter depen- 
dence (Fig. 3(c)). A possible explanation is the boundary 
condition set by the cylindrical geometry of the fiber. 
The restriction of the transverse heat flow causes heat 
accumulation and thus leads to an enlarged melt and 
evaporation area. 

The surface roughness is determined by AFM measure- 
ments at different positions with scan areas of 0.5, 2 and 
5 /im size to avoid artefacts and to minimize noise. To 
distinguish between micro roughness, which is important 
for the optical quality, and the larger scale shape which 
determines the imaging properties, the deviation from a 
polynomial fit from each scan line is used for roughness 
calculations. All measurements show comparable rough- 
ness values of a = 0.24 nm rms. Fig. 4 shows the resid- 
ual from this fit and a two dimensional power spectral 
density calculated from the data (2D PSD). Reference 
measurements on mica sheets enable us to estimate the 
noise of the measurement to be cr no [ se = 0.1 nm. Cor- 
recting for this yields a — 0.22 nm rms for the laser- 
machined structures, a value approaching that of super- 
polished optics. Scatter loss can be estimated [ ] by 
S « (47rcr/A) 2 , which predicts S ~ llppm at A = 830 nm 
for this roughness. For two-fiber cavities with state-of- 
the-art coatings, we have recently obtained a cavity fi- 
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FIG. 4: (a) AFM measurement of a 2 x 2 /im 2 area in the 
center of a laser machined area. Shown is the surface elevation 
after substracting a fitted polynomial, (b) The 2D PSD of a 
5x5 /im 2 AFM scan for the remaining height elevation. 



nesse of T = 100, 000, compatible with the surface rough- 
ness result within the uncertainty of the measured cavity 
transmission. 

While traditional mirror superpolishing methods are 
restricted to large radii of curvature in the centimeter to 
meter range [18], the method described here achieves ex- 
tremely small radii of curvature, and allows a wide range 
of geometries not limited to the range of data shown here. 
It thus opens a previously inaccessible geometry regime 
for high-quality micro-optical lenses and mirrors, start- 
ing from low-cost standard substrates. First applications 
in neutral-atom, quantum dot and molecule CQED [10- 
12, 19] have already shown the power of this method, in 
realizing microscopic, high-finesse optical cavities with 
small mode volume. Applications in cavity optomechan- 
ics and quantum information experiments with ion traps 
are under way. 
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